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On ' We continue the systematic next-to-next-to-leading order (NNLO) QCD analysis of the 

^ . CCFR data for the xF^ SF. The NNLO results for the Q^-dependence of the parameters, 

! which describe the x-dependence of this SF, are presented. 
M-i. 

D I 1. Introduction 



In the series of the subsequent papers we concentrated ourselves on the fits of 
the CCFR data for xFs SF § at the LO, NLO and NNLO level with twist-4 effects 
included through the infrared- renormalon model of Refs. . These fits were done using 
the Jacobi polynomial-Mellin moments version of the DGLAP method, developed and 
used in Refs.0- [|1^]. In the process of these fits we have taken into account the results 
of calculations of the NNLO corrections to the coefficient function of xF^ SF pT| and the 



available analytical expressions for the NNLO corrections to the anomalous dimensions 
of the nonsinglet (NS) even moments with n = 2, 4, 6, 8, 10|]l2l, supplemented with the 



given in Ref. |l[] n = 3,5,7,9 similar numbers, obtained using the smooth interpolation 
procedure, which was previously proposed in Ref. [p!3|] . 

The Jacobi polynomial method is based on the reconstruction of the SF from its mo- 
ments = x"' ^F3{x, Q'^)dx via the following equation: 

n=0 3=0 

where n = 2,3,.., Cj^\a,f3) is the combination of Euler F-functions. For simplicity we 
present here the basic formula in the case when target mass corrections and higher-twist 
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effects are neglected. More details can be found in Ref. 0], where the formulae for all 
these contributions are given. 

Making the fits with the free Jacobi polynomial parameters a,/?, we found that their 
values a ~ 0.7, /3 ~ 3 correspond to the minimum in the plane 

Using the standard model parameterization for the NS 0:^3 SF, we are now defining the 
corresponding moments at the initial scale Q^: 

MnQl)= f'x''-'A{Ql)x'^'^o)^i-a:riQl)li + ^^Ql)a;)dx (2) 







and evolve them to other transfered momentum using the solution of the renormalization 
group QCD equations for the moments, which has the following form: 

(0) 



where As{Q'^) = as(Q^)/(47r) is the MS'-scheme expansion parameter, AD{n,Q'^) = 1 + 
p{n)As{Q'^) + q{n)As{Q'^)'^ + ... comes from the expansion of the anomalous dimensions 
term of the solution of the corresponding renormalization-group equation, C^"'\Q'^) = 
1 + C'^^\n)As{Q^) + C^^\n)As{Q'^f + ... is the coefficient function of the n-th moment 
of the X-F3 SF. The numerical values of the coefficients p{n), q{n), C^^\n) and C^'^\n), 
normalized to the case of / = 4 numbers of active flavours, are given in Ref.pj^. In our 
flts we will use well-known LO, NLO and NNLO expressions of as, written down as the 
expansions in the inverse powers of /n(QVA2^^)-terms in the LO, NLO and NNLO. 

2. Presentation of the results 

Here we complete our previous analysis of the CCFR'97 data of Refs.[§-^ and concen- 
trate our attention on the extraction of the dependence of the paremeters of the model 
of Eq.(2) for xF^ SF (which, of course, are related to the parton distributions parame- 
ters). In Table 1 we present the results of our flts of the CCFR xF^ data with the cut 
> 5 GeV^ for different flxed falues of Ql with TMC taken into account, but twist-4 
contributions neglected. Several comments are in order: 

1) For the x-shape parameters the shift between LO and NLO values is larger than 
between NLO and NNLO results (exept parameter 7 ). This fact indicates the supression 
of the uncertainties of the values of the parameters as the result of taking into account 
higher order perturbative QCD effects. 

2) The numerical values of the parameters, obtained in the flts at the NNLO, are in a 



good agreement with those obtained in Ref. [0 where HT corrections were taken into 
account and the nuclear effects parameterized by deuteron model approximation of xF^*^ 



3) The effect of the increase of the values of the parameter c and the simultaneous 
decrease of the values of the parameter b, which is observed while incresing the value of 
Qo-scale, is in a qualitative agreement with the theoretical predictions of the works of 
Ref. 1161 and Ref.0. 



4) The value of the parameter 6, which governs the low-x behaviour of the xF^ SF, 
is in agreement with the recent theoretical calculation of the low-x behaviour of the NS 
structure functions ||T8|. 
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Table 1 

The Qo-dependence of the parameters of xF^^ model, extracted from the fits of CCFR'97 
data with the cut > 5 GeV^. The statistical errors are taken into account. 

Ql = b GeV^ [MeV] A b c 7 xVn.e.p. 



LO 266 ± 35 5.13 ± 0.46 0.72 ± 0.03 3.87 ± 0.05 1.42 ± 0.33 113.2/86 

NLO 341 ± 30 4.05 ± 0.36 0.65 ± 0.03 3.71 ± 0.06 1.93 ± 0.36 87.1/86 

NNLO 293 ± 29 4.25 ± 0.38 0.66 ± 0.03 3.56 ± 0.07 1.33 ± 0.33 78.4/86 
g2 = 8 GeV^ 

LO 266 ± 35 5.10 ± 0.15 0.70 ± 0.01 3.94 ± 0.04 1.23 ± 0.12 113.2/86 

NLO 340 ± 40 4.38 ± 0.17 0.66 ± 0.02 3.81 ± 0.07 1.47 ± 0.19 87.4/86 

NNLO 312 ± 33 4.42 ± 0.36 0.66 ± 0.03 3.68 ± 0.07 1.13 ± 0.31 76.5/86 
Ql = 10 GeV^ 

LO 265 ± 34 5.07 ± 0.15 0.70 ± 0.01 3.97 ± 0.04 1.16 ± 0.12 113.2/86 

NLO 340 ± 35 4.48 ± 0.15 0.66 ± 0.02 3.85 ± 0.04 1.32 ± 0.10 87.5/86 

NNLO 318 ± 33 4.50 ± 0.36 0.65 ± 0.03 3.73 ± 0.07 1.05 ± 0.31 76.3/86 
Ql = 15 GeV^ 

LO 265 ± 35 5.02 ± 0.16 0.68 ± 0.01 4.01 ± 0.04 1.04 ± 0.11 113.1/86 

NLO 339 ± 37 4.61 ± 0.42 0.65 ± 0.03 3.92 ± 0.11 1.08 ± 0.42 87.6/86 

NNLO 324 ± 34 4.60 ± 0.34 0.65 ± 0.03 3.83 ± 0.07 0.92 ± 0.29 76.6/86 
Ql = 20 GeV^ 

LO 264 ± 35 4.98 ± 0.15 0.68 ± 0.01 4.05 ± 0.04 0.96 ± 0.12 113.1/86 

NLO 339 ± 36 4.67 ± 0.16 0.65 ± 0.02 3.96 ± 0.05 0.95 ± 0.13 87.6/76 

NNLO 326 ± 35 4.70 ± 0.34 0.65 ± 0.03 3.88 ±0.08 0.80 ± 0.30 77.0/86 
Ql = 30 GeV^ 

LO 264 ± 37 4.92 ± 0.17 0.67 ± 0.02 4.09 ± 0.05 0.86 ± 0.09 113.0/86 

NLO 338 ± 37 4.72 ± 0.43 0.64 ± 0.03 4.01 ± 0.08 0.80 ± 0.35 87.5/86 

NNLO 327 ± 35 4.78 ± 0.32 0.64 ± 0.03 3.96 ± 0.08 0.67 ± 0.27 77.8/86 

Ql = 50 GeV^ 

LO 264 ± 35 4.84 ± 0.15 0.65± 0.01 4.13 ± 0.04 0.75 ± 0.12 112.9/86 

NLO 337 ± 34 4.74 ± 0.13 0.64± 0.01 4.06 ± 0.06 0.64 ± 0.14 87.5/86 

NNLO 326 ± 36 4.85 ± 0.31 0.64± 0.02 4.03 ± 0.09 0.53 ± 0.27 78.8/86 
Ql = 100 GeV^ 

LO 263 ± 36 4.73 ± 0.23 0.64± 0.02 4.19 ± 0.09 0.62 ± 0.24 112.6/86 

NLO 337 ± 37 4.73 ± 0.15 0.62± 0.02 4.12 ± 0.06 0.46 ± 0.12 87.3/86 

NNLO 325 ± 36 4.91 ± 0.28 0.63± 0.02 4.11 ± 0.10 0.36 ± 0.25 80.0/86 
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5) In the case of Qq = 5 GeV^ the NLO resuhs for the parameters b and c are in 
agreement with the outcomes of the NLO fits of the CCFR'97 data, which were made in 



Ref. |19| with the help of the standard reahzation of the DGLAP equation. Note, however. 



that the value of parameter 7, obtained in Ref.|T9|, turned out to be comparable with 
zero. It is of interest to trace the origin of this only difference with the presented in Table 
1 results. 

6) One can see that the values of A-^, which come from the LO and NLO fits, are 
rather stable to the choice of the initial scale Qq. However, the results of the NNLO 
fits are more sensitive to its variation. Indeed, despite the fact that for smaller values of 
Ql the difference between NLO and NNLO values of the parameters A, b, c, 7 are not 
large, the NNLO values of A^^^ are considerably smaller, than the NLO ones. This effect 
acreasir 

NNLO values of A.^^ are becoming rather stable. 



is decreasing while increasing Qq from 5 GeV^ up to 20 GeV"^. For > 20 GeV^ the 



MS 

7) We think that this effects is related to rather peculiar behavior of the NNLO pertur- 
bative QCD expansion of n = 2 moment. Indeed, taking into account the exact numerical 
values of the coefficients p{2), g(2), C^^\2) and C^^\2), given in Refs.gg], we find that 
the perturbative behavior of the n = 2 moment is determined by the following series: 

AD(2,Q2)c(2)(q2) = I - o,i32As{Q^) - 46.16(a,{Q^)]\ ... (4) 



where the value of the relatively large A"^ coefficient mainly comes from the NNLO term 
of the coefficient function of n = 2 moment. Thus we conclude, that it is safer to start 
the QCD evolution from the scale Qq = 20 GeV^, where the numerical value of the A'^ 
contribution in Eq.(4) is smaller, and therefore the complicated asymptotic structure of 
the perturbative expansion of n = 2 moment is not yet manifesting itself. Note, that the 
choice of this initial scale is also empirically supported by the fact that it is lying in the 
mid of the kinematic region of the CCFR data. 

8) It is interesting to note, that the nonstandard structure of the perturbative QCD 
expression for Eq.(4) is manifesting itself not only in the case of / = 4 number of flavours. 
Indeed, in the cases of / = 3 and / = 5 the expressions for the NNLO approximations of 
the corresponding solutions of the renormalization-group equations read 
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AD(2,g2)c(2)(g2) ^ 1-0.271^(Q2)- 43.51 (^^(Q2)j + .... (5) 
AD(2,g2)c(2)(g2) _ 1 + 0.126A,{Q') - A8.88( A.iQ'))' + .... 



Thus, we conclude, that the unnatural behaviour of the order 0(^4^) approximation of 
the renormalization-group improved Mellin moment cannot be avoided after changing the 
number of active flavours. 

In Table 2 we present the results of our fits of the CCFR'97 data including target mass 
corrections both without twist-4 corrections and with twist-4 contributions, fixed through 
the IRR model of Ref. as 

M„^^^ = C(n)Mr(g^)^ (6) 
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Table 2 



The results of the fits of CCFR'97 data with the cut > 5 GeV^, obtained in the case 
Ql = 20 GeV^. N^LO means the application of [0/2] expanded Fade approximants. 



[MeV] 


A 


b 


c 
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A'^[GeV^' 


xVn.e.p. 


LO 264 ± 36 


4.98 ±0.23 


0.68 ±0.02 


4.05 ±0.05 


0.96 ±0.18 




113.1/86 


433 ± 51 


4.69 ±0.13 


0.64 ±0.01 


4.03 ±0.04 


1.16 ±0.12 


-0.33 ±0.12 


83.1/86 


NLO 339 ± 35 


4.67 ±0.11 


0.65 ±0.01 


3.96 ±0.04 


0.95 ±0.09 




87.6/86 


369 ± 37 


4.62 ±0.16 


0.64 ±0.01 


3.95 ±0.05 


0.98 ±0.17 


-0.12 ±0.06 


82.3/86 


NNLO 326 ± 35 


4.70 ±0.34 


0.65 ±0.03 


3.88 ±0.08 


0.80 ±0.28 




77.0/86 


327 ±35 


4.70 ±0.34 


0.65 ±0.03 


3.88 ±0.08 


0.80 ±0.29 


-0.01 ±0.05 


76.9/86 


N^LO 335 ±37 


4.77 ±0.34 


0.65 ±0.03 


3.85 ±0.08 


0.71 ±0.28 




77.9/86 


(Fade) 340 ± 37 


4.78 ±0.34 


0.65 ±0.03 


3.85 ±0.08 


0.71 ±0.28 


-0.04 ±0.05 


77.2/86 



where C{n) = -n - A + 2/{n + 1) + A/{{n + 2) + ASi{n) {Si{n) = E"=i 1/j) and A'^ is the 
free parameter. 

The theoretical uncertainties of the NNLO results are probed by using the [0/2] Fade 
expanded approximants in Eq.(3) and using the N^LO expression for the expansion pa- 
rameter As, which depends from the calculated recently four- loop coefficient of the QCD 
/3-function pD[ . 

The problem of comparison of the results of the fits with the IRR-model estimates and 
the effect of decreasing the value of its parameter A2 at the NNLO was analysed in detail 
in Ref.|0] and Ref. [Q in particular. Here we would like to draw the attention to the fact, 
that besides the NLO value of (and thus as) is sensitive to the twist-4 contributions, 
modeled with the help of the IRR approach, the parameters of the x-dependence of the 
xFs SF (and thus parton distribution parameters) are more sensitive to the choice of Q^- 
scale, than to the twist-4 effects (see Tables 1,2 above and Table 5 in Ref. |^). Taking into 
account this effects might allow to fix theoretical uncertainties of the parton distributions 



parameters. The importance of this theoretical problem was emphasized in Ref . ||2T[ . 

To conclude, let us present the NLO and NNLO values of as{Mz), obtained in Ref. [^ 
from the results of the fits of CCFR'97 data for xF^ with twist-4 terms modeled through 
the IRR approach (see Table 2) and using the M5'-matching condition in different orders 
of perturbation theory [p^ -p^- 



NLO as{Mz) = 0.120 ± 0.003(stat) ± 0.005(syst) ± 0.004(t/ieor) (7) 
NNLO as{Mz) = 0.118 ± 0.003(stat) ± 0.005(si/st) ± 0.003(t/ieor) 

The systematical uncertainty in these results are determined by using the original CCFR 
considerations and the theoretical errors are fixed using the results of [0/2] Fade approx- 



imations fits and the proposals of Ref. |p5| , |26[l to estimate the ambiguities due to smooth 



transition to the world with / = 5 numbers of fiavours. It should be noted that the NLO 



result is in agreement with the one, obtained in Ref.|]l9[ using the DGLAF equation and 
taking into consideration the correlations of statistical and systematical uncertainties. 
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